QUASI-PERIODIC SOLUTIONS OF THE 
SCHRODINGER EQUATION WITH 
ARBITRARY ALGEBRAIC NONLINEARITIES 



W.-M. Wang 

Abstract. We present a geometric formulation of existence of time quasi-periodic so- 
lutions. As an application, we prove the existence of quasi-periodic solutions of b fre- 
quencies, b < d + 2, in arbitrary dimension d and for arbitrary non integrable algebraic 
nonlinearity p. This reflects the conservation of d momenta, energy and L 2 norm. In Id, 
we prove the existence of quasi-periodic solutions with arbitrary b and for arbitrary p, 
solving a problem that started Hamiltonian PDE. 
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1. Introduction and statement of the theorems 



We consider the nonlinear Schrodinger equation on the <i-torus T d = [0, 2Tv) d : 

d 

i—u = -Au+\u\ 2p u (p>l,peN), (1.1) 

with periodic boundary conditions: u(t, x) = u(t, x + 2mr), x G [0, 27r) d for all n G Z d . 
The integer p in (1.1) is arbitrary. The monomial nonlinearity can be replaced by 
polynomials: J2 p>po a p \u\ 2p , a p G R. For simplicity of exposition, we will only do that 
when d = 1 and po = 1 as otherwise it is integrable. 

Let u^°> be a solution to the linear equation: 

i^io) = _ Au (o). (L2) 
dt v ' 

The main purpose of this paper is to formulate two geometric conditions on u^°\ which 
are sufficient and essentially necessary for the existence of time quasi-periodic solutions 
close to to the nonlinear equation (1.1) when has finite Fourier support. These 
two conditions remain the same when u^ ' has large Fourier support and are likely to 
be relevant there as well. 

We seek quasi-periodic solutions to (1.1) with b frequencies in the form 

u(t,x) = J2 u(nj)e in ^- x , (n,j) G Z 6+d , (1.3) 

with ui G Z & to be determined. Writing in this form, a solution to (1.2) with b 
frequencies = {jl} b k=1 (jk 0) has Fourier support 

supp^°> = {(-e ifc ,j fc ),A;=l,...,6}, (1.4) 

where e Jfe is a unit vector in Z 6 . Define the bi-characteristics 

C = {(n,j) G Z b+d \ ± n ■ w (0) + j 2 = 0}. (1.5) 



C is the solution set in the form (1.3) to (1.2) and its complex conjugate in the Fourier 
space. 

Define 

C+ = {{n,j)\n ■ u {0) +j 2 = 0,j ^ 0} U {(n, 0) \n ■ u {0) =0,m< 0}, 
C~ = {(n,j)\-n- + f = 0, j± 0} U {(n, 0)\n ■ = 0, m > 0}, 
C+ n C~ = 0, C+ U C~ = C and 

C++ = {(n,j)\(n,j) = (n' - n",j'-j"), (n',j>), (n",j") G 

the set of differences and similarly C H , C h and C . 
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By supp, we will always mean the Fourier support, so from now on we write 
supp u(°) for supp etc. Let 

S = supp u {0) U supp u {0 \ (1.7) 

and 

F(u< >) = |u< >| 2 *u<°> 

the error term. For any subset S C Z b+d and / a function on T b+d of the form (1.3), 
define the restricted sum 

fs(t,x)= /("•./>'"'•"'' 

(n,j)ES 

The two conditions on are: 

(i) non intersection 

{supp F(u(°>) U supp F(«(°))} n = 0, 

(ii) non spiral 

{(n,0)|n^0}^supp {[[[\u^ \% + ] k ™ [\u^\ 2 ^{u^) 2 ] 8 ^- [l« (0) [(« (0) ) V 0) | 2(p_1) ]c-+}> 

where m > 1, /c m , fc^ > 0, s m = 0, 1; if s m = 0, then either k m or fc^ = 0. 

Remark. The cubic nonlinear Schrodinger equations fulfill (i) in any dimension d for 
all u(°> . When d = 1, it also fulfills (ii) for all u^°\ cf. sect. 5. Parameter dependent 
tangentially non resonant equations studied in [Bl, 3, EK] generically fulfill the above 
two conditions, as the geometry is "random" . 

The first condition indicates that the error term has no additional Fourier support 
on the characteristic variety C aside from and v,(°\ Since to first approximation, 
propagation is along the bi-characteristics, this is a natural condition. 

The second condition provides separation in the non convex n direction. In the j 
direction this is from convexity (cf. Lemma 2.2), fundamentally incompatible with 
translation invariance. This ensures that the lattice points on C connected by the 
convolution generated by form disjoint sets (in C) of uniformly bounded diameter. 
As we will see, this is essentially necessary. In the special case of periodic solutions 
6=1, separation is entirely provided by convexity as n is a scalar and there is positivity. 
It is useful to note that the right side of (ii) is a subset of supp {Alg\u^\ 2 } : where 
Alg\u^\ 2 = {\u^\ 2e ,£= 1,2...}. 

The main result is 

3 



Theorem 1. Let 

b 

u w (t,x) = J2 a ke lJk ' x e~ tjlt 
k=i 

be a solution to the linear Schrddinger equation (1.2) and a = {ak} b k=l G (0,5] b = 
B(0,5), jk 7^ 0, k = 1, Assume satisfies conditions (i, ii). There exist C, 

c > 0, such that for all e G (0, 1), there exists 5q > and for all 5 G (0, So) a Cantor 
set Q with 

meas {Q n B(0, 5)}/5 b > 1 - Ce c . 
If a G Q, there is a quasi-periodic solution to the nonlinear Schrddinger equation (1.1) 

b 

u (t,x) = J2 a ke tJk ' x e~ lLJkt + 0(5 2p+1 ), 
k=i 

with basic frequencies u = {u>k}\.=i satisfying 

u k =f k + 0(S 2 P +1 ), k = l,...,b,p>l. 

The remainder 0(S 2p+1 ) is in an appropriate analytic norm on T b+d . 

Remark. The above theorem also holds when there is in addition an overall phase 
m^O, corresponding to adding m to the right side of (1.1). 

The following theorems are corollaries. But since they are of a general nature, we 
state them as theorems. 

Theorem 2. There are quasi-periodic solutions of b frequencies for the nonlinear 
Schrddinger equation (1.1) for b < d + 2 in any dimension d and for any nonlinearity 
p. 

Remark. The above theorem holds on a set of a = {a k } b k=1 G (0,5] b (0 < 5 1) of 
positive measure. 

For p > 2 or d > 3, this is a new result on existence of quasi-periodic solutions. 
Basically it reflects the conservation of <i-momenta, energy and I? norm. For quasi- 
periodic solutions with cubic nonlinearity p = 1 in dimensions d = 1, 2 see [KP, Bl, 
2, GXY]. (For periodic solutions b = 1 of the wave equation, cf [BB].) 

In d > 2, (1.1) is supercritical for p sufficiently large, where there are no global 
existence results [B2]. Theorem 2 constructs a type of global solutions with precise 
control over the Fourier coefficients. 

Theorem 3. For d = 1, there are quasi-periodic solutions of b frequencies for the 
nonlinear Schrddinger equation (1.1) for arbitrary b and arbitrary nonlinearity p. 

Remark. The above theorem holds on a set of a = {cik} b k=1 G (0,5] b (0 < d <^ 1) of 
positive measure. 
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For p > 2 this is a new result. The fundamental fact in Id is that the zeroes of 
analytic functions or polynomials are isolated points, which enables us to go beyond 
integrability. 

Proof of Theorem 2. Assume b = d + 2 and supp = {(— e Jfe , jk), k = 1, ...,&}, 
b < d + 2 follow similarly. Since the right side of (ii) C supp {Alg\u^\ 2 }, it suffices 
to ensure {(n, 0)|n^ 0} supp {Alg\u^\ 2 }. 

Writing n = {n k } b k=1 leads to 

b 

J = $>fcjfc = 0, j k eZ d , (1.8) 
k=i 

which comprises of b linear equations in rik, k = 1, Since the convolution is 

restricted to C and the generator real, we have the two additional equations 



k=i 
b 

J> fc = 0. (1.10) 



k=i 

So if the determinant of the (d + 2) x (d + 2) matrix: 

/i i -i \ 

det j'i j 2 •••Jd+2 ^ 0, 

\Jl J2 '"Jd+2/ 

then n = is the unique solution to (1.8-1.10) and condition (ii) is satisfied. 
To ensure (i), we note that 

supp F(u(°>) = supp (|w (0) | 2p ) + supp w (0) , 

where the addition is element wise. Let (An, Aj) e supp (It^ ) | 2p ) and (— e m ,j m ) G 
supp w/ ). Unless (An, Aj) = (— e^/ + e m ,jV — j m ), in which case 

(-e m ,jm) + (An, Aj) = (-e k >J k >) E <S, 

the function 

/ = (-e m + An) ■ u<V + (j m + Aj) 2 = An ■ + 2j m Aj + Aj 2 

is not a constant function. It suffices to assume jk(i) are fixed for all k = 1,2..., b 
and i ^ 1. / is a polynomial of degree 2 in {jk(^)}k=i- The set <fi where / = has 
codimension 1. There are at most b 2p+1 such /. So for {jfc(l)}| =1 ^ {4>}, (i) is verified. 

Clearly similar arguments hold for 6 < d + 2 and we have proven the theorem. □ 
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A sketch of proof of Theorem 1. 

We write (1.1) in the Fourier space, it becomes 



F(u) = diag (n ■ u + j 2 )u + (u * v)* p * u = 0, (1.11) 

where (n, j) E Z b+d and v = u and to E IR 6 is to be determined. Since from now on we 
work with (1.11), for simplicity we drop the hat and write u for u and v for v. We seek 
solutions close to the linear solution of b frequencies, supp = {(— e Jfc , j k ), k = 
1, &}, with frequencies = {jl} b k=1 (jk ^ 0) and small amplitudes a = {ak} b k=1 , 
\\a\\ < 1. 

We complete (1.11) by writing the equation for the complex conjugate. So we have 

diag (n • to + j 2 )u + (u * v)* p * u = 0, 

o (1.12) 
diag (—71 • a; + j )v + (u * v)* p * v = 0, 

Denote the left side of (1.12) by F(u, v). We make a Lyapunov-Schmidt decomposition 
into the P-equations: 

F(u,v)\ zb+d \ S = 0, 

and the Q-equations: 

F(u,v)\ s = 0, 

where S is as defined in (1.7). We seek solutions such that u\s = u^°\ The P- 
equations are infinite dimensional and determine u in the complement of supp u^; 
the Q-equations are 2b dimensional and determine the frequency u = {u>k} k= i- 

This Lyapunov-Schmidt method was introduced by Craig and Wayne [CW] to con- 
struct periodic solutions for the Id wave equation. It was inspired by the multiscale 
analysis of Frohlich and Spencer [FS]. This construction was further developed by 
Bourgain to embrace the full generality of quasi-periodic solutions and in arbitrary 
dimensions d [Bl-3]. More recently, Eliasson and Kuksin [EK] developed a KAM ver- 
sion of the theory in the Schrodinger context. All the above results, however, pertain 
to parameter dependent tangentially non resonant equations. 

We use a Newton scheme to solve the P-equations, with as the initial ap- 
proximation. The major difference with [CW, Bl-3, EK] is that (1.12) is completely 
resonant and there are no parameters at this initial stage. The frequency a/ - 1 is an 
integer in Z 6 . So we need to proceed differently. 

First recall the formal scheme: the first correction 

where ^(i) ^ is the next approximation and F'(u^°\ v^) is the linearized operator 

on £ 2 (Z b+d x Z b+d ) 

F' = D + A, (1.13) 
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where 



and 



/diagKc+j 2 ) \ 

\ diag -u)+j 2 ) J v ; 

\p(u*v)* p l *v*v (p+l)(u*v)* p J ' v ; 

with ijj = uj(°\ u = and v = v^°\ 

Condition (ii) ensures invertibility of F', \\F' \\ < C(||a||- 2p ) by using bounded 
diameters of the connected sets in C and making an initial excision in a. Combined 
with condition (i), this gives 



||Au (1) || = ||Av (1) || < 0{ 



a 



\ 2 p\ 



for small ||a||. Inserting this into the Q-equations, which determine u>, we achieve 
amplitude-frequency modulation: 



HAw^H xO(||a|| 2p ) 
da 



det (^r)l x °(H a ll 2p ) >0 



ensuring transversality and moreover Diophantine on a set of a of positive measure. 
The tangentially non resonant perturbation theory in [B3] becomes available. 

The first iteration is therefore the key step and is the core of the present construc- 
tion, which is summarized in Theorem 4 at the end of section 3. We note here that 
if (ii) is violated, then there will be connected sets in C of infinite diameter due to 
translation invariance in n and typically G a(F'(u^\v^)). So both (i, ii) are essen- 
tially necessary. (Recall (i) prevents error propagation.) (ii) is an algebraic condition 
on projected loops ensuring separation in the non convex n direction in the original 
Z b+d lattice. This is a much weaker separation compared to the convex j direction. 
It suffices as we only need it for the first step. Afterwards it is taken care of by the 
modulated Diophantine frequency. 

The present theory relies on analyzing connected sets on the bi-characteristics C. For 
the cubic nonlinearity, locally this reduces to the well known resonance structure for 
the partial Birkhoff normal form [KP, Bl, 2, GXY]. So cubic nonlinearity is contained 
in the present theory as a special case. In section 5, we show this correspondence as 
well as proving Theorem 3. 

This paper is the nonlinear component of the resonant perturbation theory that we 
are developing. For the linear theory, see [Wl, 2]. 
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2. The first step in the Newton scheme— extraction of parameters 

Since we seek solutions with small amplitude a = {afc}| =1 , it is convenient to rescale: 
o-k - ► 5 1//2p ayt (0 < 5 € 1). We therefore solve instead: 

f diag (n- <jj + j 2 )u + S(u * v)* p * u = 0, ^ 
1 diag (—n ■ u> + j 2 )v + S(u * v)* p * v = 0, 

with w| supp u (o) = = a G (0, l] b = B(0, 1) and the same for v. The frequency 
lo = u(a) is to be determined by the Q-equations with the initial approximation 

We use the usual procedure of iteratively solving the P and then the Q-equations. 
The main novelty here is the first step, where u = a/ -* is a fixed integer (vector), the 
opposite of a Diophantine vector. So the usual Diophantine separation is not available. 
We overcome the difficulty by using the algebraic condition (ii) to analyze the geometry 
of the resonant structure. After the first step, with appropriate restrictions on a, 
u = u}(a) becomes Diophantine from amplitude- frequency modulation and the non 
resonant mechanism in [B3] is applicable. 

The first iteration. 

Using it(°), and a/ - 1 as the initial approximation, we first solve the P-equations. 
This requires an estimate on the inverse of the linearized operator F'(u^°\ v^) eval- 
uated at oo = u)(°>. This section is devoted to prove the following. 

Lemma 2.1. Assume u^°\ satisfy conditions (i, ii). Let a G (0, l] b = £>(0, 1) = 
B C M & \{0}. There exist C , c > 0, such that for all e G (0, 1), there exist do > and 
B e with 

meas (B e n B) < Ce c . 

If a G B\B e , then for 5 G (0, So), 

UfF'^y^)] -1 !! < 0{d~ l ) (2.2) 
and there exists (3 G (0, 1) such that 

\[FXu( \vW)}-\x,y)\<S® x - y \=e-ft lo z s U x -^ (2.3) 
for all \x — y\ > l/(3 2 . 

We need the following arithmetic fact, a proof of which is given in [B3, Lemma 
19.10], so we do not repeat it here. 

Lemma 2.2. Fix any large number B. There is a partition ofL d satisfying the 

properties 

diam ix a < B C{ \ C = C(d), 
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and 

\j-j'\ + \f-j ,2 \>B, 
if n G 7r a , n' G 7r a / u>zit/i a^a'. 

Proof of Lemma 2.1. Let P± be the projection on Z b+d onto C 1 * 1 defined in (1.6), 

on Z b+d x Z 6+d and P c the projection onto the complement. The linearized operator 
F'(u (0 \ is F' = L> + <5A with D as in (1.14) and A as in (1.15). 

From the Schur complement reduction [SI, 2], A G cr(F') fl [—1/2, 1/2] if and only 
if G a(H), where 

H = PF'P - A + PF'P C (P C F'P C - X)- 1 P C F'P (2.5) 

is the effective operator acting on the bi-characteristics C. Since || P^'P !! > 1 — 
0{8) > 1/2 and \\PF'P C \\ = 0(5), the last term in (2.5) is of order 5 2 , uniformly for 
AG [-1/2,1/2]. 

So 

H = PF'P — A + 0(5 2 ) 

in L 2 uniformly for A G [-1/2, 1/2]. To obtain (2.2), it suffices to prove || [PF'P]- 1 1| < 
0(5~ 1 ). Since PF'P = 5PAP, where A as in (1.15), let A' = PAP, it is important 
to note here that A' — A' (a) depends on a but is independent of 6. 

Let 

R > 2p diam (supp w (0) ). (2.6) 
From Lemma 2.2, there is a partition {7r a } of Z d such that 

dist (C n (Z b x tt q ), C n (Z 6 x 7r Q /)) > R (2.7) 

for a ^ a' and diam 7r a < P 1 ^ , Co = Co(d). 

Given two sets S u S 2 C Z 6+d , Si n S 2 = 0, we say that S± is connected to S^, 
Si ~ S2, if there exist si G Si, s 2 G S 2 such that A(si, s 2 ) 7^ 0. From (2.6, 2.7), 

{C n (z b x 7r a )} ^ {c n (z b x tt q /)}, 

if a 7^ a'. We say that a set £ is connected if for all s G S, s is connected to {S\s}, 
s ~ {S\s}. 

Assume S C {C fl (Z b x 7TQ,)} is a connected set and |S| > 2R c ° d . Since diam 7r a < 
R c °, there exist s, s', s 7^ s', s, s' G C + or s, s' G C _ such that s = (n,j), s' = (n',j) 
have the same Z d coordinate. So s — s' = (n — n', 0). Since S is a connected set, this 
contradicts condition (ii), unless n = n' . So 

\S\ < 2R Cod . (2.8) 
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Combining (2.6-2.8), this gives 



A' (a) = ®Aa,p{a), 

where A a ,/3 are matrices of sizes at most 2R c ° d x 2R c " d . Since A(a) is a convolution 

matrix and A(x, y) = if \x — y\ > R, there are at most Y^nLi 2 i?Ar2 < 2 R3C ° d = K 
types of A a ,(3, which we rename as Ak, 1 < k < K. 

For each fc, detAk(a) = Pfc(a) is a polynomial in a of degree at most 2 P ■ 2R c ° d < 
logK. Since for all fc, Pfc is a non constant function on B = (0, l] 6 , we have that there 
exist C, c > 0, such that for all < e < 1, 

meas {a G B\\P k \ < e, all k < K} < Ce c , (2.9) 

(cf. e.g., [Lemma 11.4, GS] and references therein). 

Since \\A k \\ < 0(1), if |detA| > e, then < 0(e~ 1 ). (The exponent is 1 

because of self-adjointness.) This proves (2.2). 

To summarize, there exists a partition of C, {S a }, \S a \ < 2R c ° d for all a, such that 
5a ^ S a >, if a 7^ a', i.e., for all s G S a , s' G 5 a /, A(s, s') = 0. Moreover 

llt^J- 1 !!^^- 1 ), 

where 

F 'sS x ^y) = F '( x ^y)^ tfx,yes a 

=0, otherwise. 
To prove (2.3), we use the resolvent expansion: 

[F , ]- 1 = [F]- 1 + [F]- 1 r[F / ]- 1 , 

where 

F = ®aF' Sa ® {n , mc diag (±n • u; (0) + j 2 ) 
V = F' — F. 

Now 

||([F] _1 r) 2m || < C(5 m ) (2.10) 

for to = 1, 2, since [F^J -1 ]?^ J -1 = 0, which is obvious when a = a' and when 
a ^ a' from (ii). (2.2, 2.10) give (2.3) for some < (3 < R~ c ° d . □ 

Using Lemma 2.1 to solve the P and then the Q-equations, we obtain the following 
result after the first iteration. Let 

AuM = - u(°), Av^ = - v(°\Au,W = w (i) _ 0,(0). 
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Proposition 2.3. Assume u^°\ satisfy conditions (i, ii). There exist C, c > 0, 
such that for all e G (0, 1), there exist So > and B e with 

meas (B e n B) < Ce c . 

If a e B\B e , then for 5 G (0, 5 ), 

\\AuU\\ = \\AvV\\ e > {p) = 0(8), (2.11) 

where p is a weight on Z b+d satisfying 

p(x) =e /3|x| , < (3 < 1 for \x\ > x 
=1, for \x\ < xq. 



\\F{u^\v^)\\ P{p) = 0{8% (2.12) 

||Aw (1) || ~0(5), (2.13) 
a ,(i) 

|det( — )|xO(5). (2.14) 

Moreover u/ 1 ) is Diophantine 

\\n-uj {1) \\ T > ^f-, n G Z b \{0}, « > 0,7 > 26+ 1, (2.15) 
|n| 7 

where \\ \\j denotes the distance to integers inW, k = K,(e) and 7 are independent of S. 
Proof. From the Newton scheme 

A(^)=[F>(0),,(0))]-^),,(0)) 
= [D' 1 +5F'~ 1 AD- 1 ]F. 

Since 

supp F(u(°\vW)r){C\S} = ®, 

from (i), 

WD-^y = o{8) 

and 

WF'-'AD-'FlU^Oil), 

where we also used (2.2). So 



ll A (^)ll^=^) 
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and 

\\F(u^,v^)y < 0(||F"||)||A \\% = 0(5 3 ). 

Using (2.3), the above two estimates hold in weighted space as well and we obtain 
(2.11, 2.12). 

From the Q equations, the frequency modulation 

A^^^F^^v^-ejJ) 
aj 

is a rational function in a of degree p + 1 for each j. So HAo/ 1 ^ 2 is controlled from 
below by a polynomial in a of degree 2(p + b). Using the same argument as in (2.9), 
we obtain (2.13). 

Similarly det(^) is a rational function in a of degree at most (p + 1)6, and is 
controlled from below by a polynomial of degree at most pb + b + 26! — 2. So (2.9) 
gives (2.14), which in turn implies the Diophantine property (2.15) for appropriate k 
as d b a = -i — d ^ ( lr-s i and the excised set in u/ 1 ) is of measure 0(5). □ 



3. The second step 

Proposition 2.3, in particular (2.15) has transformed the (tangentially) resonant per- 
turbation into non resonant perturbation by amplitude-frequency modulation. How- 
ever due to the extra power of 5, we need to proceed differently at the initial step. 
Essentially we redo Lemma 2.1 at scale N = | log 5 1 s (s > 1) with Diophantine u/ 1 ) 
replacing integer This is made possible by the fact that the resonance structure 
remains the same at this scale. 

Define the truncated linearized operator F' N {u^\ v^) evaluated at u/ 1 ) as 

F' N (uM,vM)(x,y) = F'(u^,v^)(x,y), ||x||oo < N, IMU < N 

= 0, otherwise. 

We have the analogue of Lemma 2.1. 

Lemma 3.1. Assume Proposition 2.3 holds. There exist C, c> 0, such that for all 
e G (0, 1), there exist do > and £> e with 

meas (B e n B) < Ce c , 

where B = (0, l] b . If a G B\B e , then for 6 G (0, 5 ), 

IIK^ 1 ),^ 1 ))]- 1 !!^^- 1 -) (3.2) 

and there exists (3 G (0, 1) such that 

ItFj^uWy 1 ))]- 1 ^)! < 6® x -v\ = e-P^W^ (3.3) 
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for all \x — y\> 1//3 2 . 

(3.2, 3.3) together with (2.15) put the construction in the standard non resonant 
form as put forward by Bourgain [B3], cf. also [BW]. (2.14) ensures transversality, 
in other words, continuing amplitude-frequency modulation. The extra power of 5 in 
(2.15) are compensated by the two extra powers of 8 in the remainder in (2.12). 

Proof of Lemma 3.1. Using (2.11, 2.12), it suffices to work with F' N (u^-°\ v^) (evalu- 
ated at a/ 1 - 1 ). The geometric construction in the proof of Lemma 2.1 remains valid: 

PF' N P = @ k T k {a), k < Ki(N), 

where T k are matrices of at most sizes 2R ° d x 2R ° d , R as in (2.6), 

p ( diag (n-AwW) \, X a 

Lk ~{ diag (-n ■ AwW) J + dAk ' 

and A k are convolution matrices. 

With the addition of the diagonal term, Yk is no longer a convolution matrix. 
Moreover \n\ < N = | log <5| s (s > 1) depend on 5. So we need to proceed differently 
because of uniformity considerations in estimates of type (2.9). 

Fix 

N = max(100 J R Cod , R 2C » d ). (3.4) 
For a given T^, define the support of Tk to be 

Z b+d x Z b+d D supp T fc = {(x,y)\T k (x,y) ^ 0}. 

For matrices T k: such that 

supp T k n {[-N , N ] b+d x [-N , N ] b + d } ^ 0, (3.5) 

we proceed as in the proof of Lemma 2.1. There are at most Kq (independent of S) of 
these matrices. Let P k = P k (a) = detract), we have that there exist C, c > 0, such 
that for all < e < 1, 

meas {a G B\\P k \ < 5 1+£ , all k < K } < C5 C£ . 

So Hr^ 1 !! < 0(5~ 1 ~ e ) for all k < K . It is important to note that C, c are independent 
of 5, since P k /5 are independent of 5. 

For matrices T k with k > Kq, 

supp T k n {[-N ,N ] b+d x [-N Q ,N Q ] b+d } = 

by definition. We use perturbation theory. For any F k , fix N', with \N'\ > Nq, such 
that for all (n,j) E supp T k , we can write (n,j) = (N', 0) + (n',j) with \n'\ < R Cod < 
A^"o/100. 
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Let A be an eigenvalue of 1^ with normalized eigenfunction cf>. Then 
A = (0, A0) = (0, [diag (w • Ac^ (1) ) + M fc ]0) 

= n ■ Auj{1) \<Pn,j\ 2 + S((p,A k (p) 

= N'-AJ^+ n' ■Au^\cf> n!j \ 2 + S( ( f ) ,A k( f>) 

(n',j) 

First order eigenvalue variation gives 



So 



AT' . = _/v + n( ) 



N' . no 

80,(1) > L 

AT' - 100 1 1 



Using | det(^^)| x 6 from (2.14) and taking A = 0, this gives Hr" 1 )! < O^" 1 " 6 ) 
for all K < k < Ki(N) away from a set of a of measure less than 5 e / 2 as Ki(N) < 
C(|log5| 2s ), s > 1. 

The estimates on above and the Schur reduction (2.5) give (3.2). The pointwise 
estimates follow as in the proof of (2.3), since the geometry of the resonant structure 
remains the same. Here we also used (2.11). □ 



We summarize the findings in Proposition 2.3 and Lemma 3.1 in the following 
amplitude-frequency modulation theorem. For simplicity, we use u to denote both the 
function and its Fouier series as it should be clear from the context. 



Theorem 4. Let 

b 

u w (t,x) = J2 a ke tJk ' x e~ ljlt 
k=i 

be a solution to the linear Schrddinger equation (1.2) and a = {ak} b k=1 G (0, l] b = B, 
jk 7^ 0, k = 1,...,6 and N = | log <5 1 s (s > 1). If satisfies conditions (i, ii), then 
there exist C , c > 0, such that for all e G (0, 1), there exist Sq > and B e with 

meas (B e n B) < Ce c . 

If a G B\B e , then for 5 G (0, 5 ), 

\\Au^\\ eHp) = \\Av^\\ P{p) = 0(S), 
where p is a weight on Z b+d satisfying 

p(x) =e /3|x| , < f3 < 1 for \x\> x 

=1, for \x\ < xq. 
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\\F(u^,v^)y {p) =0(8% 
||Aw (1) || xO(5), 
fit 

|det( — )|x0(o), 



ana 7 a^ 1 ) is Diophantine 



w (1) || T > j^p, n G Z b \{0}, « > 0,7 > 26+1, 

where \\ \\j denotes the distance to integers in R, k = «(e) and 7 are independent of S. 
Moreover 

MK^M^r^zotf- 1 -') (3.6) 

and 

^(u^y 1 ))]" 1 ^)! < 5 /3|a; - J/| = e -^o S 5\\x- y \ (3 _ 7) 
/or a// |x - y\ > 1//3 2 , where as in (3.1) with N = | log5| s (s > 1). 

4. Proof of Theorem 1 

Theorem 4 provides the input for the initial scale in the Newton scheme in [B3]. 
To continue the iteration, we need the analogues of (3.6, 3.7) at larger scales. This is 
achieved as follows. 

Let T = F' be the linearized operator defined as in (1.13-1.15) and the restrcited 
operator T/v = F' N as in (3.1). To increase the scale from N to a larger scale Ni, 
we pave the N± cubes with N cubes. In the j direction, this is taken care of by 
perturbation; while in the n direction by adding an additional parameter 9 G R and 
consider T(9): 

/diag(n.c + ^ + ^) \ 

\ diag (—n • uj + 2 - 9) J 

where A as defined in (1.15). 

Let iV = I logoT (s > 1) as in Lemma 3.1 and T N (9) = T N (9;u ( - 1 \v^) evaluated 
at u/ 1 ). We have the following estimates. 

Lemma 4.1. Assume u^°\ satisfy conditions (i, ii) and a G B\B e , the set defined 
in Theorem 4- Then 

\\[T N (9)]-'\\<0(S-^) (4.1) 
and there exists (3 G (0, 1) such that 



\[T N (9)]- 1 {x,y)\ < S^ x - y \ = e -^o S 5\\x- y \ 
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(4.2) 



for all \x — y\ > l/(3 2 , away from a set B^{6) C K. with 

meas B N {6) < 5 1+C£ . 

Proof. For any 9 G IR, we can always write 9 — + 59', where G Z, -| < 56' < |. 
Since ||Tjv|| < 2| log<5| 2s (s > 1), we may restrict 6 to |6| < 2\ log5\ 2s + 1. 



Then 



/ diag (n-u(°) +J 2 + ©) 
V diag (-n-wW+j 2 -e) 

/ diag (n-u + 6') \ 

+ H diag ( — n ■ Hi — 9') J ^ Ni 

where G Z 6 , G Z, w = Au/ 1 )/^ is Diophantine and is the restricted A as 
defined in (1.15, 3.1). 

Let 

/ diag (n ■ u + 6') 
V diag (-n ■ Cj - 9') 

Let P + be the projection onto the set {(n, j)\n-uj( > +j 2 +Q = 0} and P_ the projection 
onto the set {(n,j) \ —n-uj^+j 2 — Q = 0} when 6^0; when 6 = 0, use the definition 
in the proof of Lemma 2.1. Define 



H = 5 



+ A N 



P = 



P+ 
P- 



and P c the projection onto the complement as before. 

We proceed using the Schur reduction as in the proof of Lemme 2.1 and 3.1. It 
suffices to estimate PHP = ® k T k {9) as P C T N P C is invertible, \\(P C T N P C - A) _1 || < 4 
uniformly in 9 and A G [-1/4, 1/4]. 

For Tfc such that (3.5) is satisfied, we obtain 

meas {^Hr" 1 !! > S' 1 '"} < C5 ce (e > 0), 

via the estimates on detTfc as in the proof of Lemma 3.1. For other T k , we make a 
shift n ■ tb = N' ■ to + n' ■ to, where \N'\ < \ log <5 1 s , \n'\ < N , N as in (3.4). Let 
9" = 9' + N' ■ u. Then clearly 

meas {0"|dist (9",a(T k )) < 5 1+e } < 0(5 1+e ). 

Combining the above estimates on r^, we obtain (4.1). Using the geometric struc- 
ture afforded by the resonance structure, we obtain (4.2) as before. □ 

Proof of Theorem 1. This follows from Lemma 4.1 and lemme 19.13, 19.38 and 19.65 
in chapter 19 of [B3]. □ 



16 



5. Proof of Theorem 3 and the cubic nonlinearity 



Cubic nonlinearity. 

For simplicity we write u for and to for a/ -*, the solutions and frequencies of 
the linear equation. The symbols of convolution for the cubic nonlinearity are \u\ 2 , 
u 2 and u 2 . Assume (n,j) G C + ((n, j) G C~ works similarly). In order that (n,j) is 
connected to (n'j') G C, it is necessary that either 

(a) [u * v](n, j ; n'j') 7^ or 

(b) [v*v](n,j;n'j') ^ 0. 
Case (a): Since 

n ■ uj + j 2 = 0, 
n' ■ uj + j' 2 = 0, 

subtracting the two equations gives immediately 

Uk-j' k )-(j+jk) = 0, (5.1) 

where jk,jk' G Z d (/c, k' = 1, 6), are the b Fourier components of u. 
Case (b): Since 

n ■ uj + j 2 = 0, 
-n'-u + j' 2 = 0, 
adding the two equations gives immediately 

(j-3k)-(j-jk') = 0, (5.2) 

where jk,jk' £ ^ d (k, k' = 1, 6), are the b Fourier components of u. 

(5.1, 5.2) are precisely the well known resonant set for the partial Birkhoff normal 
form transform in [KP, Bl, 2, GSY]. (5.1, 5.2) describe rectangular type of geometry, 
(i) is satisfied for the cubic nonlinear Schrodinger in any d. In Id, (5.1, 5.2) reduce to 
a finite set of 2b lattice points in Z: j = ±j k , k = 1, b. The set {=Fefc, ±jk} is the 
only connected set in C and (ii) is verified for all u^°K 

The following proof extends the existence of quasi-periodic solutions in Id of any 
number of frequencies b to arbitrary algebraic nonlinearities p. It shows that it is not 
necessary to know explicitly the locations of the resonances, which is the case for cubic 
nonlinearity, except that they are isolated. 

Proof of Theorem 3. 
Let 

r+ = supp [\u^] = {(An,Aj)}, 
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where 

An = - ^2pkk'(ek - e k '), 
Aj = ^Pkk'Uk ~ jk'), 
Pkk> > 0, ^Pkk> < P, k, k! = 1, b, 
{(-e fc ,ifc)}Li = supp n (0) , 

and 

r- = supp [l^l 2 ^- 1 ^^) 2 ] = {(An, Aj)}, 

where 

An = - ^2pkk'(e k -e k >) + (e K + e K >), 
Aj = ^Pkk'Uk - jk') ~ Uk+3k'), 
Pkk' > 0, ^Pkk* <p-l,k, k', k, k' = 1, b 

Assume {j k } b k=1 * s sucn that (An, 0) ^ r + U T~ for n^O, which is always possible 
from (5.3, 5.4). Assume (n,j) G C + ((n,j) G C~ works similarly). Then in order for 
(n,j) to be connected to (n'j') G C via (An, Aj), at least one of the following two 
equations has to be satisfied as in the cubic case: 

2jAj + Aj 2 + An ■ u = 0, if (An, Aj) G T + (5.5) 

or 

2j 2 + 2jAj + Aj 2 - An-u = 0, if (An, Aj) G T". (5.6) 

Let {j,j + Aj} be the Z coordinates of the set of connected pairs determined by 
(5.5, 5.6). Let D = { j - j', j + Aj - j', j - j' - Aj', j + Aj - j' - Aj', Aj ^ Aj'} be 
the set of differences. The constant functions in D are only generated by pairs where 
Aj = j m ±j n (and similarly for Aj'), in which case (j, j + Aj) = (Tjnjm)- These 
are the terms present in the cubic nonlinearity. 

For the non constant functions /, / = implies P = for some polynomial P in 
{j k } b k= i- Let 4> be the zero set of P. It has codimension 1. There are at most 8b 2p 
such polynomials. Choose {j k } b k=1 ^ {</>}• Condition (ii) is then verified. Fulfilling (i) 
as in the proof of Theorem 2 concludes the present proof. □ 

Remark. Theorem 3 is coherent with the result in [B4] proving existence of almost 
periodic solutions with subexponential decay amplitudes for the parameter dependent 
non resonant Id quintic nonlinear Schrodinger equation. 



(5.3) 



(5.4) 
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